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ABSTRACT 

Pulsar wind nebulae are a prominent class of very high energy (E > 0.1 TeV) Galactic sources. 
Their 7-ray spectra are interpreted as due to inverse Compton scattering of ultrarelativistic electrons 
on the ambient photons, whereas the X-ray spectra are due to synchrotron emission. We investigate 
the relation between the 7- and- X-ray emission and the pulsars' spin-down luminosity and charac- 
teristic age. We find that the distance-independent 7- to X-ray flux ratio of the nebulae is inversely 
proportional to the spin-down luminosity, (oc -E~ 1,9 ), while it appears proportional to the character- 
istic age, (oc t 2,2 ), of the parent pulsar. We interpret these results as due to the evolution of the 
electron energy distribution and the nebular dynamics, supporting the idea of so-called relic pulsar 
wind nebulae. These empirical relations provide a new tool to classify unidentified diffuse 7-ray sources 
and to estimate the spin-down luminosity and characteristic age of rotation powered pulsars with no 
detected pulsation from the X- and 7— ray properties of the associated pulsar wind nebulae. We apply 
these relations to predict the spin-down luminosity and characteristic age of four (so far unpulsing) 
candidate pulsars associated to wind nebulae. 

Subject headings: pulsars : general — radiation mechanisms: non-thermal — supernova remnants — 
X-rays : stars — gamma rays: observations 



1. INTRODUCTION 

Pulsar Wind Nebulae (PWNe) arise when the wind 
ejected from a rotation powered pulsar is confined by 
the pressure of the surrounding medium, be it their 
supernova remnant or compressed interstellar gas (see 
iGaensler fc Sla nc 2006, for a review). The Galactic sur- 
vey performed by the H. E. S. S. exp eriment {High Energy 
Stereoscopic System, iHintonl l2004f) has detected several 
PWNe making them a prominent class of Very High En- 
ergy Galactic so urces ([Aharonian et al.1 l2006eb I Gallant! 
l2007t [Funk 2003). In addition to the classical investiga- 
tions through radio and X-ray astronomies, Very High 
Energy 7-rays (VHE 7-ray s, E > 0.1 TeV) provide a 
new probe of the physical con ditions in PWNe (e.g., 
Ide Jaeer fc Diannati-Atall2008[ ). 

The PWN broad-band radiation is most likely due to 
electron-positron pairs of the pulsar wind generated close 
to the magnetosphere. The wind flow is ultrarelativis- 
tic (bulk Lorentz factor Tw ~ 10 6 in the Crab Neb- 
ula; [Kennir&iCoroniH [1984110), until it experiences a 
strong shock, where electrons are accelerated. After the 
shock, the flow speed is sub-relativistic at the outer edge 
of the PWN. Depending on the radiation mechanisms 
at work, the electrons can produce photons in different 
energy ranges: while synchrotron radiation yields pho- 
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tons with energies up to several MeV, inverse Compton 
scattering of the ambient photon field can produce high 
energy photons, up to tens of TeV. 

The electrons responsible for the PWNe 7-ray emission 
(here after 7-ray electrons) are likely less energetic than 
those generating the X-ray one (X-ray electrons), their 
synchrotron radiation lying at infrared, optical, or ul- 
traviolet frequencies. For typical nebular magnetic field 
intensities (B ~ 1-100 (J.G), synchrotron photons with 
energy ~1 keV are produced by electrons with Lorentz 
factor ~0.3-3 xlO 8 . The Cosmic Background Radia- 
tion, the dust-scattered light, and the starlight provide 
the target photons for inverse Compton scattering, with 
typical photon energies around 10 -3 eV, 10 -2 eV, and 1 
eV, respectively. In the Thomson regime, photons with 
energy ~1 TeV are produced by electrons with Lorentz 
factor ~0.1-3 xlO 7 . Due to their different energies, the 
cooling time of the X-ray electrons is smaller than the 
one of the 7-ray electrons. Therefore, the X-ray emis- 
sion traces the recent history of the nebula, whereas the 
7-ray emission traces a longer history, possibly up to the 
pulsar birth. The different lifetime of the electrons, to- 
gether with the interaction with the ambient medium, 
can lead to the significant projected angular separation, 
sometimes measured between the peaks of the 7- and X- 
ray br ightness profiles (e.g. G18.0-0.7. lAharonian et ail 
I2006dh . Since the source of injected electrons, the pulsar 
rotational energy loss rate dubbed spin-down luminosity, 
decreases as time goes by, we expect a different evolution 
of the 7- and X-ray luminosities, following the particle 
aging and the pulsar spin-down. 

In this paper we address first the correlations between 
the PWN VHE 7-ray luminosities (1-30 TeV) and their 
X-ray luminosities (2-10 keV) with the spin-down lu- 
minosities, E, and the characteristic ages, r c , of their 
pulsars. Next we consider the behaviour of the ratio be- 
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tween the gamma and X-ray luminosity as a function 
of the pulsar spin-down power and age. These relations 
are discussed in the frame of an evolving electron energy 
population. 

2. OBSERVED CORRELATIONS 

In Table [1] we report a sample of the identified 
PWNe observed by the H.E.S.S. experiment. We fur- 
ther included six candidate PWNe, selecting unidenti- 
fied H.E.S.S. diffuse sources located near young and en- 
ergetic pulsars, with r c < 100 kyr and E > 10 35 erg 
s _1 . These parameters are defined as E = 4ir 2 IP/P 3 
and t c = P/2P, where P is the pulsar spin period, 
P its derivative, and / = 10 45 gm cm 2 the moment of 
inertia. We calculated E and t c using the P and P val- 
ues reported in the Australia Telescope National Faci l- 
ity (ATNF) pulsar catalogue 7 (jManchester et al.l l2005h . 
The 7-ray fluxes, Fy, are derived from literature and 
computed in the 1-30 TeV energy band, with statistical 
errors estimated with standard Montecarlo propagation 
technique. The lower energy value corresponds to the 
highest observed detection threshold. The upper value 
of 30 TeV reduces the bias of possible unmeasured high- 
energy cut-offs. The unabsorbed X-ray fluxes, Fx, have 
been derived from literature based on X-ray imaging ob- 
servatories, and converted in the 2-10 keV energy band. 
The lower energy is chosen in order to minimize the con- 
tamination by possible thermal components due to the 
pulsar or supernova remnant. When it was possible to 
single out the PWN from the pulsar component, only the 
PWN flux is reported. 

We investigated the relations between the different lu- 
minosities and the pulsar parameters, using the data col- 
lected in Table [TJ The 7-ray luminosities, L 7 , do not ap- 
pear correlated with the pulsar spin-down luminosities E, 
nor they do with the characteristic ages r c , as shown in 
Fig. [TJ(top panels). This is at variance with the observed 
PWNe X-ray luminosities, for which a scaling relation is 
apparent with both E and r c (Fig. [TJ middle panels). 
The weighted least square fit on the whole dataset yields 

log 10 L x = (33.8 ± 0.04) + (1.87 ± 0.04) log 10 E 37 . (1) 

All the uncertainties are at la level, and E = E37 x 
10 37 erg s . The Lx — E scaling is known for 
the pulsars as well as for the PWNe. This scal- 
ing was firstly noted by ISeward fe Wand (|1988h . fur- 
ther, iBecker fc Triimpen (|1997h investigate a sample of 
27 pulsars with ROSAT, yielding the simple scaling 
£ j y (o.i-2.4keV) — 10 -3 -E. A re-analysis was performed 
by iPossenti et alj (|2002t ). who studied a sample of 39 
pulsars observed by several X-ray observatories, account- 
ing for the statistical and systematic errors. They found 
log 10 L x = (-14.36 ± 0.01) + (1.34 ± 0.03) log 10 E, a re- 
lation harder than Eq. ([TJ) . However, they could not 
separate the PWN from the pulsar contribution. A 
better comparison can be done with the results from 
iKargaltsev fe Pavlov! (J200S[) , who recently used high- 
resolution Chandra data in order to decouple the PWN 
and the pulsar fluxes. Indeed, taking E, r c , and Lpwn 

7 http: //www. atnf . csiro . au/research/pulsar/psrcat 



in the 0.5-8 keV energy band 8 from their Tables 1 and 
2, we obtained as fitted values log 10 £x(o.5-8keV) = 
(34.02 ± 0.05) + (1.46 ± 0.04) log 10 E37 for their whole 
sample, and log 10 Lx(o.5-8kcV) = (34.26±0.03) + (1.87 ± 
0.01) log 10 E37 restricting the fit only to the sources also 
present in our sample. The latter is compatible in the 
terms of slope with Eq. ([1]), and the slight difference in 
normalization can be due to the different energy band. 
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Fig. 1. — 7-ray luminosity, X-ray luminosity, and 7- to X-ray 
flux ratio versus pulsar spin-down luminosity, E (left column), 
and characteristic age, r c (right column). Filled and open cir- 
cles stand for identified and candidate PWNe, r espectively. The 
upper-limit for the flux ra tio of PSR B1706-44 (Aharonia n et al.l 
l2005al ;lR omani et al. 2005) is reported with an arrow. Also shown 
are the best-fit curves for identified PWNe (dotted lines), and for 
the whole sample (dashed lines). 

X-ray sources of our whole dataset also show a depen- 
dence of Lx on r c , with a best-fit relation 

log 10 L x = (33.7 ± 0.04) - (2.49 ± 0.06) log 10 r 4 , (2) 

where r c is in units of ye ars. The Lx — t c scal - 
ing was already noted by IBecker &; Triimperl (|1997l ) 
and IPossenti et al.l (|2002f h Also in this case we com- 
pared our fit to the one derived using the whole 
IKargaltsev fc Pavlov! (|2008h dataset, which results in 



logio L 



10 ^X(0.5-8keV) 



(34.29±0.01) - (2.03±0.01) log 10 r 4 



for their whole sample, and log 10 L 



10 ^X(0.5-8keV) 



(34.23± 



8 The X-ray luminosity reported in Kargaltsev & Pavlov ( 2008) 
for Kes 75 wa s corre cted according to the distance measured by 
ILeahv fc Tianl pool) . 



Evolution of the 7- and X-ray luminosities of PWNe 
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TABLE 1 

Properties of Wind Nebulae observed with H.E.S.S. and associated Pulsars 



Source 


Associated 


F° (1-30 TeV) 
10 -12 erg cm -2 s _1 


F b x (2-10 keV) 


T c 


E 


Distance 


References 


Name 


Pulsar 


erg cm -2 s — 1 


kyr 


erg s _1 


kpc 





v^raD 




-IV 


oU 


w 


(lb) 


in 
Z. 1U 


X 


1U 


8 


1 1 


•1 ft 

i.O 


X 


m38 

1U 


1 nq+0.11 

^•^-o.ii 


1 Q 
l,z,o 


Vole 


dcd Rns'^'^ 

r ori, DUooo- 


— <±o 


70 


(10,1 


(IbJ 




X 


1 n- 


11 


1 1 

1 1 




X 


1U 


n OQ7+0.019 




K3 in Kookaburra 


PSR J1420- 


-6048 


14.5 


(1.6) 


(2.9) 


1.3 


X 


iir 


12 


13 


1.0 


X 


10 37 


5.6+g-g 


7,8,9 


MSH 15-52 


PSR B1509- 


-58 


20.3 


(1.1) 


(4.1) 


2.86 


X 


lo- 


11 


1.6 


1.8 


X 


10 37 


5 2+ 14 


10,11,12 


G18.0-0.7 


PSR B1823- 


-13 


61 


(7) 


(12) 


4.4 


X 


ur 


13 


21 


2.8 


X 


10 36 


q+0.4 

j.y_ .4 


13,14,9 


G21. 5-0.9 


PSR J1833- 


-1034 


2.1 


(1.1) 


(0.5) 


4.0 


X 


10" 


11 


1.9 


3.4 


X 


10 37 


n n+0.4 

o -°-0.5 


15,16,9 


AX J1838.0-0655 


PSR J1838- 


-0655 


18.0 


(2.7) 


(3.6) 


1.0 


X 


10" 


12 


23 


5.5 


X 


10 36 


6-6±°;| 


17,18,19 


Kes 75 


PSR J1846- 


-0258 


2.3 


(0.6) 


(0.5) 


2.27 


X 


10" 


11 


0.73 


8.1 


X 


10 36 




15,20,21 


HESS J1303-63lt 


PSR J1301- 


-6305 


12 


(1.2) 


(2.4) 


6.2 


X 


10" 


14 


11 


1.7 


X 


10 36 


a c+1.2 
6 - 6 -l.l 


22,23,9 


HESS J1616~508t 


PSR J1617- 


-5055 


21 


(3) 


(4) 


4.2 


X 


10- 


12 


8.1 


1.6 


X 


10 37 


6 7+ - 7 

4 7+0-5 
4 -'-0.5 


17,24,9 


HESS J1702-420 1 " 


PSR J1702- 


-4128 


9.1 


(3.4) 


(1.8) 


6.0 


X 


10" 


15 


55 


3.4 


X 


1Q 35 


17,25,9 


HESS J1718-385t 


PSR J1718- 


-3825 


4.3 


(1.3) 


(0.9) 


1.4 


X 


10" 


13 


90 


1.3 


X 


10 36 


3 6 +0 - 4 


26,27,9 


HESS J1804-216 1 " 


PSR B1800- 


-21 


11.8 


(1.6) 


(2.4) 


6.8 


X 


10" 


14 


16 


2.2 


X 


10 36 


0+0.4 

3-S_ . 5 


17,28,9 


HESS J1809-193 1 " 


PSR J1809- 


-1917 


19 


(4) 


(4) 


2.3 


X 


10" 


13 


51 


1.8 


X 


10 36 


c+0.4 
J -°-0.5 


26,29,9 



■(Candidate sources. a 7-ray fluxes, statistical, and systematical errors. When not stated in the original papers, the systematic errors were assumed 
at the typical value of 20% as in I Aharonian et al.l d2006eD. b Er rors are conservatively estimated at 20%. 

References.- (1) Aharonian ct al. 2006c; (2) Willingalc et al. 2001; (3) Trimble, V. 1973; (4) Aharonian ct al. 2006a; (5) Manzali ct al. 2007; (6) 
Dodson ct al. 2003; (7) Aharonian ct al. 2006b; (8) Ne et al. 2005: (9) Manchester ct al. 2005; (10) Aharonian ct al. 2005c; (11) Gacnslcr ct al. 
2002; (12) Gacnslcr et al. 1999; (13) Aharonian et al. 2006d; (14) Gacnslcr ct al. 2003; (15) Djannati-Atai et al. 2007; (16) Slanc et al. 2000; 
(17) Aharonian ct al. 2006c; (18) Gotthclf & Halpcrn 2008; (19) Davics ct al. 2008 (20) Holland ct al. 2003; (21) Leahy & Tian 2008; (22) 
lAharoniarTe t al. 2005d; (23) XMM public da ta ar chive; (24) iKargaltsev et all 120081 ; (25) IChang et al.ll2tl08l (26) lAha ronian ct af] 120071 ; (27) 
Hinton et al. 2007: (28) Kargaltscv et al. 2007; (29) Kargaltscv & Pavlov 2007. 



0.02) - (2.60 ±0.02) log 10 r 4 restricting the fit only to the 
sources also present in our sample. 

The lower panels of Fig. [T] refer to the 7- to X-ray 
flux ratio F 1 jFx- There is a clear anticorrelation be- 
tween F~ ( /Fx and E, spanning over four decades in E 
and seven decades in F 1 jF x (Fig. [TJ bottom left panel). 
Considering only the identified PWNe, the correlation 
coefficient is r = —0.7 ± 0.2; including also the candi- 
date sources, the anticorrelation is more significant, with 
r = —0.84 ± 0.09. The best-fit including only the identi- 
fied sources yields 

log 10 F 7 /F x = (0.47±0.05)-(1.87±0.07)log 10 £ 3 7. (3) 

For all the data points, it results 

log 10 F 7 /F x = (0.57 ± 0.04) - (1.88 ± 0.05) log 10 E 37 , 

compatible within the errors with the relation obtained 
using only the identified sources. 

The 7- to X-ray flux ratio is also found to correlate with 
the characteristic age r c (Fig. [TJ bottom right panel), 
with a correlation coefficient r = 0.7 ± 0.2 for identified 
PWNe only, and r — 0.75 ± 0.f3 for the whole sample. 
The ordinary weighted least square fit only for the iden- 
tified PWNe yields 

log w F 7 /F x = (0.70±0.06) + (2.21±0.09)log 10 r 4 , (5) 

and for all the data points 

log 10 F 7 /F x = (0.89 ±0.04) + (2.14 ± 0.07) log 10 r 4 . (6) 

One should note that these correlations are based on 8 
identified sources, and are consistent when the 6 candi- 
date sources are considered. More 7-ray detections may 
improve their significance. 

3. DISCUSSION 



We found the 7- to X- ray luminosity ratio L^/ L x = 
F 1 /F x to be anticorrelated with the spin-down lumi- 
nosity E and correlated with the characteristic age r c . 
Formally, such dependencies are driven by the scaling 
law of the X-ray luminosity L x , which increases with E 
and decreases with t c , since the values of L 1 were found 
uncorrelated with the pulsar parameters. However, the 
F 1 jF x is a distant-independent indicator which relates 
two electron populations, differing by energy and age. An 
evolution in the PWN broad-band spectrum is pointed 
out by Eq. ([5]), which implies L 7 > L x after ~5 kyr 
from pulsar birth. Therefore, the 7-ray emission remains 
efficient around L 7 ~ 10 33 -10 35 erg s _1 , while the X-ray 
luminosity decreases by a factor ~10 6 in 10 5 yr following 
the pulsar spin-down. 

Such a broad-band spectral evolution can be ex- 
pected on the basis of the PWNe lep tonic model 
(|Kennel fc Coroniti|[l984bl : IChevalieH[200a ). In a PWN, 
the source of the injected electrons is the pulsar spin- 
down luminosity, E. The total injection rate of the elec- 
trons can be written: 

E 



N= r 27TT V ( ? ) 

1 w rn e & (1 + ct) 

where the magnetization parameter <j sets the fraction of 
the spin-down luminosity converted in kinetic energy of 
the wind. The whole spin-down luminosity is converted 
in p article kinetic energy for a <C 1, as for the Crab Neb- 
ula (|Kennel &: CoronTtJ 1984a. b). For sake of simplicity, 
we assume a constant wind Lorentz factor Tyy upstream 
the s hock. E decreases in time as (e.g.. iPacini fc Salvatil 
[19731 ) 

^ = 7 T-^F' (8) 

(1 + t/t dec y 

where Eq ~ 10 38 -10 40 erg s _1 is the spin-down luminos- 
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TABLE 2 

PWNe hosting a neutron star without detected pulsations 



Source 


F 7 (1-30 TeV) 


F x (2-10 keV) 




E* 


Name 


erg cm~ 2 s _1 


erg cm -2 s — 1 


kyr 


erg s 


G313.+0.1 Rabbit 


1.0 X lO" 11 


7.3 X 10~ 12 


~6 


~ 1.5 x 10 37 


GO.9+0.1 


3.3 x 10~ 12 


5.8 x 10~ 12 


~4 


~ 2 x 10 37 


G12.8-0.0t 


1.3 x lO" 11 


9.2 x lO" 12 


~6 


~ 1.5 x 10 37 


HESS J1640-465+ 


9.3 x 10" 12 


1.0 x lO" 12 


~13 


~ 5 x 10 36 



t Candidate sources. * Predicted va lues. References.- HESS J 1418/G313.+0.1 Rabbit : 
l Aharonian eFall j2006bTl . INg et all pX)5ft ; HESS J1747 -281/G0.9+0.1: lAharonian et a.1.1 <l2Q05bl1. 
IPorquet et al.l [120031 ); HESS J1813- 178/G12.8-0.0: lAharonian et al.l d2006d1. IHelfand et al.l 
(2007); HESS J1640-465/G338.3-0.0: lAharonian ~a l. ( 2006c) ." iFunk et a l. (2007). 



ity at the pulsar birth, tdec ~ 100-1000 yr is a character- 
istic decay time, t is the time elapsed since pulsar birth 
(to = 0), and j3 = (n+ l)/(n— 1), where n is the braking 
index. In the following, we assume a pure dipolar mag- 
netic field torque, i.e. n = 3. As the braking indices in- 
ferred from the measurement of the per iod and its deriva- 
tives are significantly smaller than 3 (|T7ivingston e et all 
2007), we dealt with a generic n (see App. [SJ, and found 
that the results derived from Eq. ([5]) are unaffected by 
the choice of n. 

Since it depends on E, also the particle injection rate 
N decreases in time. Therefore, the total number of par- 
ticles 



N(x I E{t')dt' = E t de 
'o 



t 



t + tde 



(9) 



reaches a constant value N oc E^tdec for t 3> tdec, and 
the particle supply by the pulsar becomes negligible. 

The electron energy distribution n{E, t) accounting 
for particle injection and radiati ve losses evolves accord- 
ing to the kinetic equation (e.g.. | Ginzburg fc SvrovatskiH 
[l96l IBlumenthal fc Gould II1970D : 



(10) 



where Q = Q(E,t) is the particle distribution injected 
per unit time, and P = P{E, t) is the radiated power per 
particle with energy E. The normalization of n(E, t) is 
set by N via the injection rate: N(t) — j Q(E,t)dE. 

At energies for which the radiative losses are negligible, 
the number of particles n(E,t) with energy E at time t 
has the same profile of the injected distribution Q(E) 
with a normalization set by N. Therefore, 



n u (E,t) oc / E{t')dt' = E Q tde 
Jo 



t + tde 



(11) 



where u stands for uncooled. As in Eq. |9]), a constant 
value n u (E, t) oc E$ tdec is reached for t 3> td ec - 

The effect of the radiative losses is to limit the accu- 
mulation of particles at a given energy. After an energy- 
dependent cooling time t c (E), the particles with initial 
energ y E have radiat ed a significant fraction of their en- 
ergy (|Chevalie"rll2000[ ). Accounting for pitch- angle aver- 
aged synchrotron and inverse Compton in the Thomson 
regime energy losses, the cooling time can be written as 



t c (E) 



4(1+0 e 4 J E B 2 



24.5(l+£r 1 7f 1 5 5 - 2 



kyr, 
(12) 



where — E /(m e c 2 ) is the particle Lorentz factor, and 
£ = Uph/Us, with Uph and Ub the photon field and mag- 
netic field energy densities, respectively (7^ = 77 x 10 7 , 
B = B5 X 10~ 5 G). When the photon field is provided 
by the Cosmic Background Radiation (U p h = 0.26 eV 
cm~ 3 ), the synchrotron radiation is the main cooling 
process (£ < 1) if B > 3 /iG. This condition is gener- 
ally fulfilled in PWNe as the equipartition magnetic field 
intensity ranges in B ~ 1 — 100 /iG. 9 Eq. ([12"]) shows that 
the cooling time of 7-ray radiating particles, t CJ1 is one 
order of magnitude longer than that of the X-ray radiat- 
ing particles, t c x, e.g., for B = 10 /iG, t C7 ~ 8-250 kyr, 
and t c x ~ 0.8-8 kyr. By comparing t cl and t c x with 
the average characteristic ages of pulsars in TeV PWNe, 
the 7-radiation is produced by long-lived electrons trac- 
ing the time-integrated evolution of the nebula, even up 
to the pulsar birth, whereas the X-ray emission is gener- 
ated by younger electrons, injected in the last thousands 
of years. 

Only the particles injected since the last t c (E) years 
will contribute to n(E,t). Eq. (fTTj) is accordingly modi- 
fied: 



n c (E, t) oc 



E{t')dt' = - 

-tc \ l 



E t d ec f c 



where c stands for cooled. 
E t dec t c t' 2 for t > max(i c , 
Et c using Eq. ©. 



tc + tdec) {t + tdec) ' 

This implies n c (E,t) oc 
tdec), and hence n c (E,t) oc 



4. CONCLUSIONS 

Eqs. (TlT|) and (fTB"]) describe the time evolution of a 
particle populations in two regimes, uncooled and cooled. 
Such an evolution is exemplified in Fig. [2] for the pop- 
ulations of particles producing 7-rays, n 7 , and X-rays, 
nx- After the initial rise, both the particle populations 
reach a plateau (t > tdec)- The decrease begins when the 
evolution time is greater than the cooling time. As in 
general t c x < tc-y, the X-ray emission fades long before 
the 7-ray one. 

As the characteristic ages of the pulsars powering 
a VHE 7-ray PWN are in the range 1-20 kyr, likely 
t c x < t c < t CJ . Accordingly, the population of the X-ray 
electrons, nx, is likely to be in the cooling regime, i.e., 
it decreases. The scaling laws nx oc t~ 2 and nx oc E 
of Eq. (fi"3"|) support the trend observed in the data, see 



9 In radiation-dominated environment, like the Galactic Cen- 
ter, the inverse Compton can contribute to the cooling. In 
this case, the Klein-Nis hina regime should be taken into account 
UManolakou et al.ll2007l) . 
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Fig. 2. — Time evolution of the number of particles radiating 
in VHE 7-rays, n 7 , and in X-rays, nx (solid lines), and of their 
ratio (dashed line). Pulsar birth is at t = 0. Initial conditions for 
the pulsar spin-down luminosity are E = 10 39 erg s" 1 and t dec = 
100 yr. Both curves are normalized to their maximum value. After 
the initial rise, both particle populations reach a plateau. The fall 
begins at t greater than the cooling time, which is assumed to be: 
t c x = 2.6 kyr for X-rays, t cl = 25 kyr for 7-rays (for a magnetic 
field intensity B = 10 fiG, and a Lorentz factor of 7-ray radiating 
electrons 7 = 10 7 ). 

Eq. (fTJ). At variance, the population of the 7-ray elec- 
trons, n 7 , is in the uncooled regime, the asymptotic limit 
of Eq. (fTTj); this explains the lack of correlation between 
7-ray luminosity L 1 and E. Finally, Eqs. ((11)) and (| 13|) 
for t c x < t c < t cl imply a ratio n 7 /nx oc t 2 oc E^ 1 . 
Since the luminosities are roughly proportional to the 
population of radiating particles, finally one gets 

L^/Lx oc t 2 oc E^ 1 , 

to compare with the best-fit empirical relations L 7 /Lx oc 
t 2 - 2 and Lj/L x oc E~ lM , see Eqs. © and ©. Though 
the outlined model does not correctly predict the slopes, 
not surprisingly in being simplified, it highlights the con- 
current roles of the evolving pulsar injection and of the 
radiative losses in producing the observed trends. 

The scattering around the relations for F 1 /Fx reflects 
the lack of correlation between L 1 and E. Environmen- 
tal factors can affect the 7-ray luminosities, like the lo- 
cal energy density of the ambient photon field, or the 
interaction with the surrounding medium causing an en- 
hancement in the magnetic field. Also, unmeasured pul- 
sar properties such as the magnetic field, its orientation 
with respect to the spin axis, and the initial spin period 
might affect the pulsar wind properties. 

We stress that the relations presented here are derived 
under several assumptions, the most important of which 
being that X-ray and 7-ray emitting particles are in dif- 
ferent cooling regimes, cooled for X-rays and uncooled 
for 7-rays. However, the Lorentz factors ranges of the 
two populations get closer, and they can even overlap, if 
the nebular magnetic field is very high, on the order of 
£?=170 /iG. On the other hand, in the case of a young 
nebula with a very low magnetic field, the X-ray elec- 
trons may not have reached the cooling regime, leaving 
the 7-ray production to the low-energy freshly injected 



electrons. Hence, PWNe with a very weak magnetic field, 
like 3C 58 (|Slane et al.ll2008l ). or po ssibly with a unusu- 
ally st rong one, as reported lately bv lArzoumanian et al.l 
(2008) for DA 495, could represent outliers to our derived 
relations. These regimes can be properly taken into ac- 
count through numerical modelling of the kinetic equa- 
tion (Eq. [10j). Another important assumption is a uni- 
form and constant magnetic field: indeed high resolution 
imaging observations of several PWNe show a dynamical 
and structured nebular morphology (Gaens ier &: Slanel 
2006). The evolution of the average magnetic field is 
complicated by the interaction with the supernova ejecta, 
which is expected to occur after a few thousands of years 
since pulsar birth, causing global oscillations of th e mag- 
netic field intensity (e.g.. iBucciantini et al.l [20031. One 
should note that the cooling time is not well defined if it is 
comparable to or longer than the time scale of variation 
of the magnetic field. The cases of patchy or evolving 
magnetic field are further sources of scattering around 
our relations. 

Given the limitations discussed above, the empirical 
relations in Eqs. ([3]) and ([5]) provide a new tool to esti- 
mate the spin-down luminosity and characteristic age of 
a pulsar lacking detected pulsation from the 7- and X-ray 
properties of the associated PWN. For the four candidate 
pulsars in Table 2, we used F~/Fx to predict E and r c . 
Taking into account the average scattering (average ab- 
solute ratio) around the best fit relations, Eqs. Q and 
(O, one should expect an uncertainty of a factor ^2.5 
for E and ^2.3 for r c considering only the eight identi- 
fied sources. On the other hand, considering Eqs. ([4]) 
and (|6|), and including also the candidate sources, the 
uncertainties are ~2.2 for E and ~2.6 for r c . 

The correlations for F 1 /Fx hold also after including 
the candidate sources. The pulsars possibly associated 
to the candidate sources are mostly older Vela-like pul- 
sars, with 8 x 10 3 yr < t c < 9 x 10 4 yr, and 3.4 x 10 35 erg 
s _1 < E < 1.6 x 10 37 erg s _1 . Due to the pulsar ages, the 
electrons had the time to be advected far from the pulsar, 
producing the offset between the 7-ray emission centroid 
and the pulsar po sition, the process leading to t he so- 
called relic PWNe (jdOageE&^Iann^rSMlQQl) . The 
fact that all the confirmed associations contain younger 
pulsars is hence not surprising, as the positional coinci- 
dence is one of the main identification criteria. If the 
identification of candidate sources with relic PWNe is 
confirmed, the 7-ray luminosity would persist up to 10 5 
yr, with remarkable time-integrated energy channeled in 
radiation (~ 3 x 10 45 -3 x 10 47 erg). 
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APPENDIX 

PARTICLE POPULATION INJECTED BY A PULSAR WITH GENERIC BRAKING INDEX 
By adopting Eq. © for a generic braking index n > 1, Eqs. (jTTJ) and (fT5)) are so modified: 



n u (E,t) oc / E(t')dt' = 



i 



tde 



t + t de 



p-i 



and 



n c (E,t) oc 



E(t')dt' = 

t-t a P 



1 



1 - 



tr 



^dec ~t~ t 



1-/3 



- f 



(Al) 



(A2) 



For t > t dec Eq. (|A1|) yields n u oc E t dec /((3 - 1), while for t 3> max(i c , i dec ) Eq. (|A2|) yields n c oc E{t)t c . As in 
the case of the dipolar magnetic braking, in the uncooled regime most of the radiating particles has been injected in 
the early phases, whereas in the cooled regime the particle population is limited by the cooling time and follows more 
closely the current spin-down rate. 
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